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Abstract. Coherent states in the time-energy plane provide a natural basis 
to study adiabatic scattering. We relate the (diagonal) matrix elements of the 
scattering matrix in this basis with the frozen on-shell scattering data. We 
describe an exactly solvable model, and show that the error in the frozen data 
cannot be estimated by the Wigner time delay alone. We introduce the notion 
of energy shift, a conjugate of Wigner time delay, and show that for incoming 
state p{Hq) the energy shift determines the outgoing state. 



1. Introduction 

Scattering from a slowly changing scatterer is described, to leading order, by a 
time independent scatterer frozen at the scattering time . Although this seems 
like stating the obvious, it turns out that in trying to make precise how accurate 
this approximation is, one encounters both conceptual and technical difficulties. 
Our aim is to describe these difficulties and explain how they are resolved. 

One conceptual difficulty is to understand what the frozen S matrix — a function 
of energy and scattering time — means. Strictly speaking, a function of both time 
and energy is in conflict with the uncertainty principle. A wave that is sharp in 
energy will have an ill-defined scattering time and conversely, a wave with a well- 
defined scattering time is ill-defined in energy. What, then, is the meaning of the 
frozen S matrix? 

The resolution of this problem is related to the fact that the adiabatic limit 
naturally leads to different parameterizations of time, and the right parameteriza- 
tion has small uncertainty. Specifically, the physical time t will parameterize the 
intrinsic "fast" dynamics and has the usual time-energy uncertainty h. The slow 
variation in the external conditions will be parameterized by s. We refer to the 
latter as epoch. Since the epoch often plays a role of a parameter it is convenient to 
choose s dimensionless. The two parameterizations are related by s = ut, with lj 
a slow frequency — the adiabaticity parameter. The epoch-energy uncertainty then 
takes the form Ss Se ^ fko and so arbitrarily small in the adiabatic limit. 

Coherent states provide a convenient basis to analyze the semi-classical limit 
12, ^. Semi-classical analysis is traditionally about the h ^ limit, but is equally 



valid when h is fixed (and henceforth set equal to one) and lo 0. Here we 
introduce coherent states labelled by points in the time-energy plane, with time 
being the scattering time. As we shall see, the frozen S matrix approximates the 
diagonal matrix elements of the dynamical scattering matrix in such coherent states. 
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This reconciles the time-energy uncertainty with the frozen scattering data. In a 
further step, matrix elements of the frozen S matrix can be approximated by the 
on-shell data. 

Another thorny issue that we address is when a description in terms of frozen 
data is meaningful and how accurate it is. The question can also be rephrased as 
a question about the intrinsic time scale relevant to scattering. If r denotes this 
time scale, then lot is the error in the frozen data and lot <C 1 characterizes the 
adiabatic regime. 

The Wigner time delay t^{E, s) conveys information about the time the particle 
spends near the scatterer. It is a function of the energy E and scattering epoch 
s. It is tempting to hope that r might be estimated by Tw{E,s), but there is no 
compelling argument for doing so. One cannot argue on the basis of dimensional 
analysis alone, since t^ and y^r^, with dot a derivative with respect to the epoch 
and prime with respect to the energy, give additional and independent time scales. 
In fact, since Wigner time delay is a comparison of the arrival time at a faraway 
point, relative to the time of arrival in the free dynamics, it is not even positive- 
definite. This suggests that it cannot quite capture r, which is more closely related 
to the "dwell time" near the scatterer. 

The way to determine t is to consider the error in approximating the scattering 
data by the frozen data. The error is, to leading order, proportional to the adia- 
baticity parameter uj. Since the error is, in general, complex, we identify r with the 
absolute value of the error divided by lo. Calculating the error, to leading order in 
w, is no harder, and reminiscent of, calculating the scattering in the lowest order 
of the Born approximation. 

We shall see that, to leading order, the adiabatic time scale r can be estimated 
from the scattering data and the derivative of the Hamiltonian H with respect to 
the epoch, Eq. ( [7.3[ ) below, but not from the Wigner time delay alone. We show 
this by considering an exactly soluble model where the dynamical S matrix can be 
computed explicitly. 

We introduce the energy shift operator E. This is a measure of the energy change 
in time dependent scattering and is a natural dual of the Wigner time delay. As 
we shall see, in the case that the incoming state is p{Ho), the outgoing state is 
p{Ho — Lo£). In the adiabatic limit, the energy shift can be approximated by the 
frozen energy shift, which is related to the logarithmic derivative of the on-shell 
scattering matrix with respect to the epoch, Eq. ( |4.2| ). The energy shift then gives 
a handle on the exchange of energy ^ and the pumping of charge in adiabatic 
scattering ||]. 



2. Elements of scattering theory 

Scattering theory is a comparison of dynamics: One is the actual dynamics 
generated by the time dependent H{t) = Hg, {s = ojt), the other is a fiducial 
dynamics generated by a time independent Hamiltonian Hq- The Hamiltonian Hq 
is the generator of dynamics for which there is trivial scattering and the S matrix 
is the identity. 

The results of this section are true in general, without taking the adiabatic limit 
UJ ^ 0. We shall assume that H and Hq admit good scattering. Namely, we 
assume the existence of wave operators and the unitarity of the S matrix. For 
exphcit conditions on Hq and H(t) that guarantee this see e.g. 15 1. 
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2.1. The wave operator. Let U{t",t') and Uo{t",t') = Uo{t" ~ t') denote the 
evolution from time t' to i", generated by H{t) and the time- independent Hq re- 
spectively. 

Definition 2.1. The wave operators, based at epoch s, are defined by the (strong) 
limit 

(2.1) n±{s;H,Ho)= hm U{t,t')U^{t' ~ t), (s^Lut). 

t' — >±CXD 

The existence of the limit, and the equation of motion imply 

Proposition 2.2. The dependence of the wave operator on the base point s satisfy 
the differential equation 

(2.2) -iu;n±{s)^H,n±{s,H,Ho)-n±is,H,Ho)Ho. 

As we shall presently see, the notion of wave operator based at epoch s is only 
interesting in the case of a time dependent H{t). 

2.2. The frozen wave operators. The frozen Hamiltonian Hg is time indepen- 
dent so U{t",t') = e*^=^(* ^* \ in this case and Q±{sq, Hg, Hq) is independent of 
the base point sq = ujIq. This follows from the existence of the limit in Eq. (2.1) 
since t' — > ±oo is the same as t' ~ to —^ ±oo. To stress this we write D,±{Hs, Hq). 
From Eq. (2^) then follows the standard intertwining relation of time-independent 
scattering theory: 

Corollary 2.3. The wave operators Q±{Hs, Hq) relating the frozen Hamiltonian at 
epoch s and Hq are independent of the base point, and intertwine the two dynamics: 

(2.3) Hgn±iHg,HQ)^n±iHg,HQ)HQ. 

2.3. The dynamical S matrix. The (dynamical) scattering matrix based at 
epoch s is defined by 

(2.4) Sd{s;H,HQ) ^nl{s;H,HQ)n-is;H,HQ). 

The S matrices based on different points in time are all related by conjugation 
generated by the free evolution. Namely: 

Proposition 2.4. Suppose that the wave operators exist. Then 

(2.5) S4s;H,HQ) = e-'""'S40:H,Ho)e'""', (s^cot). 

This follows from [/(s, t)n±{s; H, Hq) = fi±(s; H, HQ)e-'"«^''-*l Under a change 
of the reference Hamiltonian, say to the frozen Hamiltonian Hg, 

(2.6) Sdis;H,HQ) = nl{Hg, Hq)S4s; H, Hg)n^{Hg, Hq). 

2.4. The frozen S matrix. In the frozen S data the epoch is decoupled from 
time. As such it can also be studied using time independent methods, which are 
normally quite powerful [p^ . Its basic properties are in marked contrast with that 
of the dynamical S matrix, namely: 

Corollary 2.5. The frozen S matrix 

(2.7) Sf{Hg,HQ) ^nl{Hg,HQ)n_{Hs,HQ) 

is independent of the base point. It depends on the freezing time parametrically 
through Hg. 
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2.5. The on-shell S matrix. Hq provides a basis that spans the Hilbert space of 
scattering states. Let \E,j) denote the generahzed eigenvectors of Hq: 

(2.8) Ho\E,j)=E\E,j), iE,j\E',j')^S{E-E')6,,j>. 

E is the energy and j labels the scattering channels. Sf commutes with Hq, by 
Eq. (U), hence 

(2.9) iE,j\Sf{H,,Ho)\E\j')^SiE^E')S,,-is,E). 

Sjj'{s,E) is the on-shell scattering matrix. Note that in the frozen Hamiltonian 
the physical time is decoupled from the epoch, which now has been relegated to 
the role of a parameter. The on-shell scattering matrix therefore is not in conflict 
with the uncertainty principle. 

3. The energy shift 



By taking the s-derivatives of Eq. (2.5) one gets 
iioSd{s)Sd{s)^ =Ho-Sd{s)HoSd{s)^ = [ilo, 5d(s)]5^(s) 

(3.1) = [Ho,Sdis)-SfiHs,Ho)]Slis). 

This equation may interpreted as follows. If we think of Hq as the asymptotic 
observable associated with the outgoing energy, then i?o,in — Sd{s) Hq Sdis^ rep- 
resents the asymptotic observable Q corresponding to the incoming energy. This 
motivates calling 

(3.2) Sdis) = iSdis)Sl{s). 

the operator of energy shift. 

The energy shift vanishes for time independent scattering, as it must. It gives a 
handle on changes in (certain) quantum states. By the functional calculus applied 



to Eq. (3.1), for any function p: 

(3.3) Sdis)p{HQ)Sd{sy = p{Ho-Lo£sis)). 

This is interpreted as follows: If p{Ho) is the incoming state, then the corresponding 
outgoing state is p[Hq — uj£s{s)Y The energy shift is a first order quantity in the 
adiabaticity parameter and, as we shall see, it can be approximated, to leading 
order by the frozen data. This then gives a handle on the outgoing state p to first 
order in the adiabaticity parameter. 

Proposition 3.1. The energy shift based on time s is conjugate to the energy shift 
based on time zero 

(3.4) £dis)^e'"'>'£d{0)e-'""\ (s^cut). 



This follows directly from Eq. (2.5) and Eq. (3.1 



4. The problem of adiabatic scattering 

The dynamical S matrix has qualitatively different properties from the frozen S 
matrix: The dynamical S matrix has no freezing time — It does not "know" when the 
incoming wave is going to hit the scatterer. It does depend however, by conjugation, 
on a choice of a base point. In contrast, the frozen S matrix is independent of the 
choice of a base point and depends non-trivially on the freezing time — The frozen 
scattering data for one epoch know nothing a-priori about the corresponding data 
at any other epoch. 
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Matrix elements of the scattering matrix carry information about the time that 

the wave is near the scattcrcr. For such matrix elements, the adiabatic limit can be 
expressed in terms of the corresponding frozen matrix elements. However, the intro- 
duction of wave packets promotes the epoch from playing the role of a parameter, 
to that of real, albeit slow, time. One then needs to confront the uncertainty prin- 
ciple. We do that by considering matrix elements between coherent states labelled 
by points in the energy time plane. 

4.1. The Wigner time delay. The Wigner time delay is defined in terms of 
the on-shell scattering matrix. When this definition is transcribed to the frozen, 
on-shell, S matrix it reads 

(4.1) T^{s,E) = -iS'{s,E)S''{s,E). 

Prime denotes partial derivative with respect to the energy. With this definition, 
the Wigner time-delay is a Hermitian matrix. 

4.2. The frozen energy shift. For the frozen, on-shell, Hamiltonian one can 
associate a matrix of energy shift which is a natural conjugate of the Wigner time 
delay: 

(4.2) 8{s,E) = iS{s,E)S^{s,E), 
where dot denotes derivative with respect to the epoch. 

4.3. Time scales. The frozen on-shell S matrix defines several time scales. Among 

them: and the (dimensionless) time scale £~^. The coherent states provide us 
with yet another time scale related to the time-width of the coherent states. One 
of the problems of adiabatic scattering is to study the relation between these time 
scales and the time scale r such that wr <C 1 characterizes the adiabatic regime. 

5. Time-Energy Coherent states 

5.1. The role of dispersion. For a particle moving on the line, its energy and 
the time that it crosses the origin are canonical coordinates. One can therefore 
construct energy-time coherent states in analogy with the usual phase space coher- 
ent states. The explicit construction, however, depends on the dispersion law. For 
linear dispersion the construction is particularly simple. 

Consider a classical particle with dispersion law e{p) moving freely on the line. 
The velocity of the particle is e'{p) so the time of passage through the origin is 
t = ^7^. Time-energy are (local) canonical coordinates since 

(5.1) de A dt = dq A dp. 

The global aspects of the energy-time phase space can be complicated. For example, 
for a free (massive) particle, with quadratic dispersion e{p) = p^ the energy-time 
phase space is made of two copies of the half plane e > depending on the direction 
of crossing of the origin. 

A simpler situation is obtained in the case of linear dispersion, e{p) = p. There 
is now no ambiguity in the direction of crossing and the energy-time phase space is 
again the plane. The map {q,p) ^ (e, t) is, in fact, the identity 

(5.2) e = p, t = -q. 

The usual coherent states are then also the coherent states on the energy-time 
plane. 
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5.2. Coherent states for linear dispersion. The time-energy coherent states 
are 

(5.3) |t,e;£)=e^(*™)|ge), [P,X] = -j, 
with Gaussian: 

1 _ 

(5.4) {p\ge) = ^^e ^. 



They have the following properties |12[| : 
A The states |i, e;e) are normalized. 

B |f, e; e) have Gaussian localization in time and energy near the point 
{t, e) with width 

(Se ~ £, 6t ^ 6s ^ —. 

£ e 

Hence to plays the role of h in the epoch-energy plane. 
C Ho is the generator of shifts of the coherent states: 

e-*^°*>, e; £) = e-'*'^/2 \t-t',e;e). 

D The overlap of coherent states is: 

{t,e,e\t ,e ,e) ^ e e ^ e * 2 . 

E The coherent states give a resolution of the identity 

dt de 



2n 



|t,e;e) {t,e;s\ = 1. 



F The scalar product between coherent states and the eigenstates of Hq 
P is 



6. Scattering between channels with linear dispersion 

Linear dispersion approximates the low energy physics of electrons in one dimen- 
sional channels provided the Fermi energy is large. The price one pays is that the 
"ultraviolet" properties are pathological. In particular, the spectrum is unbounded 
below and this then leads to certain anomalies which must be correctly interpreted. 
With linear dispersion one can also solve certain models with interacting electrons 

In the following we shall study adiabatic scattering for non interacting particles 
with linear dispersion. The particles move on a collection of lines and are allowed 
to "hop" from one line to the other and scatter. Each line serves as an incoming 
and outgoing channel since the flow on it is uni-directional. An example with two 
channels is shown in Fig. |l|. (Such models bear some resemblance to Schrodinger 
operators on graphs The Hilbert space is ®"^]^L^(]R), a finite direct sum. j 
labels the scattering channels. Hq is then 

(Hotp) {x,j) = -iij'{x,j), xeR, l<j<n. 

For the interaction one may take, for example, 



His) - Hq)^p] {x, j) = ^ Vjj, (x, s)i;{x, j') 
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with Vjji hermitian, and compactly supported. Alternatively, one may consider 
finite rank perturbations. 




Figure 1. A network of two channels. Each channel is chiral 
and lets particles propagate to and from infinity, according to the 
arrows. The circle denotes the region where the channels are cou- 
pled. 

6.1. A soluble model. Here we describe a simple, time dependent, model for 
which the calculation of both the dynamical and frozen scattering matrices is re- 
duced to quadrature. 

Consider scattering on the line with 

Ho = P^-iV, Hs = P + fis)V, s = ujt 

with (V4'){x) = v{x)tp{x) a potential (multiplication operator) which is sufficiently 
regular and short range so that J \v{x)\ dx, J \xv{x)\ dx < oo . The model has 
one channel and should not be confused with the 2-channel example pictured in 
Figure |l|. 

To calculate the dynamical S matrix note that 

(6.1) n{t,t') ■.= u{t,t')Uo{t' -t), 

satisfies the Volterra type equation: 

(6.2) ^^}Ml^,f^^t')nit,t')vit~t'), nit',t')^i, 

with V{t) the (backward) free Heisenberg evolution of the potential, i.e. 

(6.3) Vit) := Uoit) V Uoi~t). 

Since Hq — P is the generator of shifts, Vit) is the shifted potential: 

(6.4) {V{t)^) (x) = vix - t)^ix). 

In particular, Vit) at different times commute, and the solution of the Volterra 
type problem is given simply by 

(6.5) nit, t') = e-' -fo fi'^-^t-mndt" ^ ^^^^^ 
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From the definition of the wave operators based on time s, Eq. (|]^), we obtain for 
the dynamical wave operators: 

From this we obtain for the dynamical scattering matrix 

(6.6) Sdis,H,Ho) = e-^-f— 

The dynamical scattering matrix, as well as the wave operators, are local gauge 
transformations, i.e. multiplication by a function of position, of modulus one. 

The wave operators and the S matrix reduce to the frozen ones upon replacing 
the function /(s — cut') by its frozen value /(s), hence: 

(6.7) Sf{Hs,Ho) = e-'^(^) -f-^ ^^'^ = e'^f^^^ -f-^ ^+(*) 

where 2V+{t) — V{t) + V{—t). Sf is just a number, not a function of position. 

The frozen scattering matrix provide very little information on the potential 
v{x), for it depends one just one number — the total weig ht of the potentia]|]. The 
dynamical S matrix, in contrast, provides independent information about the po- 
tential for each value of s. 

Since the frozen S matrix is independent of the incident energy, the Wigner time 
delay vanishes identically in this model: t^, — 0. The (frozen) energy shift is just a 
real number (a multiple of the identity) 



^/ /(*) / v{x)dx. 

J —OO 

In contrast, the dynamical energy shift, is the multiphcation operator: 
(6.8) £d= / f{s-ujt')V{t')dt'. 



6.2. The on-shell scattering matrix and coherent states. For later purposes 
we shall need the matrix elements of the frozen S matrix. Since Sf commutes with 
Hq, the matrix elements are independent of t and are related to the on-shell matrix 

by 

{t,e,j;e\Sf{H,,Ho)\t,e,f;e) = J dE 3,^. {s, E)e-^^ 

(6.9) - Sjj>is,e) + Oie^dEES). 

The estimate is obtained by observing that Sjji (s, E)~Sjji (s, e) does not contribute 
to the integral to first order in E — e. Since 



{OeeS) 



_ 2 I • ' 



(with prime denoting the derivative with respect to the energy) we see that the on- 
shell S matrix approximates the diagonal entries of the frozen S matrix, provided 
the Wigner time delay and its energy dependence are both small: 

(6.10) e2(r^ + |r:|)«l. 



^This is in sharp contrast with scattering problems where Hq is the Laplacian 
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7. The adiabatic time scale r 

In this section we compute, to leading order, the time scale r relevant to adiabatic 
scattering. This time scale is defined so that wr ^ 1 characterizes the adiabatic 
regime in the sense that the frozen scattering data approximate the dynamical 
scattering data. 

There are two results in this section, one positive and one negative. The positive 
result says that, at least to leading order, r can be computed from time independent 
quantities alone, Eq. (7.3) below. The negative result is that r cannot be computed 
from the on-shell scattering matrix and its derivatives. In particular, the Wigner 
time delay alone does not determine r. 

Using Eqs. (2T5|,2.6,p77) and property 5.2.C one finds 

(7.1) {t, e, j; e\ (5^(0; H, Ho) - Sf{Hs,Ho)) \t, e, £> - 

(0, e, j; £| nl{H,,Ho){Sd{s; H, H,) - l)l]_(i/„ i/o) |0, e,i; e) . 

The correction to the leading order of the S matrix can be approximated by an 
analog of the Born series |jl^ : 



/oo 
-OO 



dt'. 



Since Hg^^^ti—Hg is supported near the origin, only small t' contribute to the matrix 
elements in Eq. (7.1). More precisely, this depends only on the time localization 
property of either the bra or the ket. We can therefore approximate Hg+^v — -ff s ~ 
ujt'Hs- Using property 5.2.C 

e-'"^'n-{Hs,Ho) |0, e, j; e) = e-'^*/^ ^_(^Hs,Ho) \t, e, j; e) 

we finally get 

(7.2) {t, e, r, e\ (5^(0; H, Ho) - Sf{Hs,Ho)) \t, e, /; e) « -zc^r(e, s; e) 
where 

/•oc 

(7.3) r(e,s;£)= / {t' ,e,j;e\nl{Hs,Ho)Hsn^{Hs,Ho)\t\e,j;e) t' dt' . 



T(e, s;e) involves the frozen wave operators and the rate of change of the Hamil- 
tonian at the epoch s. In particular, one can use methods of time-independent 
scattering theory to compute it. It is in general complex. The adiabatic time scale, 
T — \T{e,s;e)\ is a measure of the error, cut <C 1 then clearly characterizes the 
adiabatic regime. 

Propagation estimates can, and have been, used | |T^ to bound the error in the 
frozen data. These estimates yield bounds on r. 

7.1. Example: the soluble model. For the case of one channel scattering with 
H{s) + fis)V, by Eq. i^^^^^) 

S4s;H,Ho)-Sf{Hs,Ho) = ( ^-^ {fis-^t)-ns))v(t) dt _ A Sf{Hs,Ho) 



(7.4) « iu;f{s)[ tV{t)dt] Sf{Hs,Ho) 
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The adiabatic time scale r is, in analogy with Eq. (]7.2D, the multiplication operator 



(7.5)r«-/(s) 



tVit)dt\ SfiH,,Ho) = -fis) 
2V^it)^V{t)-V{-t). 



tV^it)dt\ SfiHs,Ho 



By Eq. (S.7) the frozen S matrix only depends on V+, while the error only depends 
on V-. Since V- and V+ are independent this shows that the error term in the 
adiabatic expansion cannot be estimated in terms on the frozen scattering data 
alone. 

Combining Eqs. (3.9) and (7.2) we obtain a relation between matrix elements of 
the dynamical S matrix and the on-shell S matrix: 

(7.6) 

{t,e,j;e\Sd{0;H,Ho) \t,e,f;e) = Sjj,{s,e) + 0{e\Tl + |r;|) + wr(e, s; e)) , 

8. The energy shift 

The energy shift is a first order quantity, nevertheless, it is determined, to leading 
order, by the frozen data: 



(8.1) 



(t,e,j;e|£rf(0)|t,e,/;e) « z (5(s, e)5t(s, e))^.^.,, (s ^ cot) 



We first remark that (7.2 7.3) generalize to off-diagonal matrix elements, i.e., the 
time t in the ket \t,e,j';e) may be shifted to t + At (resp. t' + At in (7.2)) while 
leaving the bra unchanged. By the translation property of coherent states, property 
5.2.C, multiplication by Hq can be traded for derivative with respect to time. Hence 

{t,e,j;e\ [Ha,Sd{0; H, Hq)] \t + At,e,j';e) 

= idt {t,e,j;e\Sd{Q;H,Ho)\t + At,e,j';e) 

« tdt {t,e,j;e\SfiHs,Ho)\t + At,e,j';e) 

(8.2) = iuj{t,e,j;e\Sf{H,,Ho)\t + At,e,j';e). 

In principle, the order of the error in the frozen data in the passage from the second 
to the third line does not determine the order of the error in derivatives, but this 
can be justified in the present case. The last identity in the equation above can be 
seen from 

{t,e,j;e\Sf{Hs,Ho) \t + At,e,j';e) 



dES, 



' — 7^ 



g-j(At)e/2g-j(At)iS^ 



We then multiply (|]^) with the complex conjugate of the mentioned generalization 
of (7.2) and integrate over At using property 5.2.E. The result then heuristically 
follows from Eq. (3T) and the statement for £f analogous to (|6.S|) .p] 

The energy shift plays a role in the theory of adiabatic quantum pumps. In 
particular, the pumped charge, the entropy production and noise generation in 
quantum pumps can all be expressed in terms of the energy shift |^ . It is remarkable 
that basic properties of adiabatic quantum pumps can be understood, to leading 
order, in terms of the frozen scattering data alone. 



^ An alternate derivation of i.l can be made, more directly, starting from the rhs of Eq. 
and using Born's expansion. 



3.1 
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